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Abstract 



This paper reports thermopower and conductivity measurements through 
the metal-insulator transition for 2-dimensional electron gases in high mobil- 
ity Si-MOSFETs. At low temperatures both thermopower and conductivity 
show critical behaviour as a function of electron density. When approaching 
the critical density from the metallic side the diffusion thermopower appears 
to diverge and the conductivity vanishes. On the insulating side the ther- 
mopower shows an upturn with decreasing temperature. These features have 

much in common with those expected for an Anderson transition. 
PACS: 71.30.+h, 73.40.-c 
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I. INTRODUCTION 



Scaling theory of non-interacting, disordered, electron gases predicts that no metal- 
insulator transition (MIT) occurs in 2 dimensionsEHi as temperature T ^ 0. Nevertheless, 
what appears to be an MIT has been observed (at finite, though low T), first in n— Si- 
MOSFETsS and more recently in many other 2-dimensional (2D) hole and electron gases.i 
In the particular case of Si-MOSFETs, the transition is most clearly visible in high-mobility 
samples, roughly /i > Im^/Vs. As the electron density, tt., is varied, there is a particular 
value, no, above or below which the resistivity p shows metallic or insulating temperature 
dependence respectively. For the present purposes we will use as a working definition that 
negative dp/dT indicates an 'insulator', and positive dp/dT at the lowest temperatures we 
can reach corresponds to a 'metal' (possible deviations from this definition and the conse- 
quences will be mentioned later). At n > uq and not too close to no, metallic behaviour 
is visible over a wide range of T, roughly T < 0.5Tp where Tp is the Fermi temperature. 
The decrease of p in the metallic state for high mobility Si-MOSFETs samples is typically 
a few orders of magnitude larger than can be accounted for by electron-phonon scattering. 
There is no accepted explanation of this behaviour as yet. On the insulating side p is found 
to increase exponentially as T is reduced, typically showing good agreement with variable 
range hopping models. 

Most previous work on these systems has focused on p, though measurements of the 
compressibilityi have also appeared recently. The present paper presents experimental data 
on thermopower, 5*, and conductivity, a = 1/p, at temperatures down to 0.3 K. When we 
extract the diffusion thermopower, S'^, we find that it diverges when approaching uq. At no 
there is an abrupt change in behaviour of the thermopower with lower densities showing an 
upturn in as T is decreased. We also attribute this to S'^. 

In addition we find that a exhibits critical behaviour around no which appears to be 
largely consistent with the behaviour of S'^. Earlier, a scaling behaviour was describedS for 
the temperature dependence of p{T) over a temperature range ~ (0.05 — 0.3)Tp. In contrast, 
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the present experiments are concerned with the hmiting, low temperature behaviour of 
a. As aheady mentioned, p decreases with temperature when n > but it saturates 
at low temperatures. (We are unable to say what happens below 0.3 K and return to 
this point later). Our low temperature saturation values of a, say ctq, show a power-law 
critical behaviour as a function of {n/riQ — 1). Similar critical behaviour of a around Uq 
has previously been reported in 2D p-type GaAs0 and lower mobility Si-MOSFETs.i The 
observed characteristics of S"^ and do are strongly reminiscent of those expected for an 
Anderson MIT in but such a transition should not occur in 2D. 

We have reported some of our preliminary results, particularly on the thermopower, in 
conference form.El The present paper extends the thermopwer results to a wider range of 
temperatures and densities, and extends the data into regions of high sample resistivities 
which were inaccessible to us earlier. It is interesting to note that, unlike resistivity, S is 
measured with no current flow in the system and hence no self heating occur. We also 
present the results of systematic studies on p on a number of samples; these were primarily 
taken to augment our understanding of the thermopower, but they are also interesting in 
their own right. 



II. SAMPLE AND EXPERIMENTS 

The main sample used for the present p and S measurements (Sample 1) is the same 
as that described in a previous paperlli and the general techniques used to measure S can 
also be found there. This sample has no = 1.01 x 10^^ m^^ (as defined as above) and a 
peak mobility p = 1.82m^/Vs at T = 0.3 K. S and p have been measured as a function of 
T, down to about 0.3 K, at many different values of n. We have also analyzed independent 
p(T, n) data for two other samples over the same range of T, Sample 2 from the same wafer 
with uq = 0.96 X 10^^ m^^ and with peak p = 1.96m^/Vs (again at T = 0.3 K) and Sample 
30 with no = 0.96 x 10^^ m ^ and peak p = 3.6m^/Vs. All connections to the 2D gas had 
isolation resistance > 50 Gfl and all leads into the cryostat were well shielded and filtered 
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against rf interference. For the sample leads, commercial, bulkhead, rf rejection filters were 
used covering a range 1 MHz - 10 GHz, and these were augmented by extra, series inductors 
to extend the range down to about 50 kHz. 

Measurements were made at dc. With the sample in the metallic state [n > no), a 
Keithley 182 digital voltmeter usually gave the best compromise of input bias current, input 
impedance and noise. As previously,0 we obtained spurious voltages from the sample even 
with no temperature gradient or current present. We traced this to the input bias current 
of the voltmeter. With any amplifier, the input bias current combined with the source 
impedance produces spurious (offset) voltages which must be separated from the true signals. 
Spurious voltages have been noticed previously by many people and probably have a number 
of causes, but we verified the consistency of the explanation in the present case by measuring 
the 2-terminal resistance of the sample over a wide range of conditions, as well as measuring 
the input bias current of the amplifier. The bias current of the input amplifier of our 
Keithley 182 was 15 pA for sources up to about IMQ but it increased to 30 pA at higher 
source impedances. 

When the system was metallic the offset voltages were always rather small and not 
usually too T dependent, especially at lower temperatures where the thermoelectric signals 
are also small, so they were not a serious problem. For example, at n = 2.1 x lO^^m^^, the 
sample is metallic and has a resistance of ~ 20 kQ at low temperatures. This gives an offset 
voltage of ~ 30 nV. Offsets were taken into account by taking the thermoelectric signal to 
be the difference between the measured voltages with the temperature gradient present and 
not present, ensuring that the mean sample temperature, and thus resistance, was the same 
in both cases. 

On the other hand, when the sample was in the insulating state its resistance was much 
higher and strongly T— dependent. In this case the offset voltage rapidly overwhelmed the 
thermoelectric signal as the temperature was lowered. To reduce this problem, an amplifier 
with input bias current < 1 pA and input impedance > lO^^Q was used. With some averaging 
it had a resolution of 0.1 for source impedances of less than a few hundred kf2 rising to 
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about 1 /iV at 10-20 MQ,. However, when used with the sample there was always an apparent 
offset current of ~ 4pA independent of sample voltage, i.e., offset signals were produced 
corresponding to a dc current through the sample of this magnitude. We were unable to 
determine the origin of this current. It was unaffected by changing the rf filters in the leads, 
or by rerouting some of the wiring, particularly that carrying the gate voltage. When the 
offset signals become large, one would like to have some check that the method of extracting 
the signal used above was reliable. Fortunately, it was found that if the sample was grounded 
at different contacts, the offset current always flowed in a direction towards ground. Thus by 
measuring the thermoelectric signals with the source or drain grounded in turn, the offset 
voltage reversed sign but the thermoelectric voltage did not. This gave two independent 
measurements and is the origin of the error bars on our data at low temperatures with the 
sample in the insulating state (see Fig. (4) later). The width of the error bars typically 
increases very rapidly as T is lowered because of the strong variation of p. 

Finally we provide some information on the temperature difference AT across the sample 
at various temperatures so that the signal voltages AV = SAT may be estimated under 
various conditions. The electric field in the sample may also be estimated given the length 
was 2.5 mm. Between 0.7 and 7K, AT ^ 100 mK. Below 0.7 K AT dropped steadily until 
by 0.4 K it was about 30 mK and at 0.34 K about 15 mK. 

III. RESULTS AND DISCUSSION 

In the metallic region n, as obtained from the Hall constant, was found to be a linear 
function of gate voltage Vg, and it is believed to follow approximately the same dependence in 
the insulating region Jll at least close to uq. On different cool downs from room temperature 
there were small shifts in n at fixed Vg of up to 0.04 x 10^^ m~^, but if the sample was held 
at 77 K between cooldowns, as was done for most of these measurements, the density and 
other properties were very reproducible. In general the absolute uncertainty in n should be 
no worse than ±0.1 x 10^^ m^^ in the region of Uq and 2% at high n; the reason for the 
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larger uncertainty near no is the fact that the Hall resistivity shows non-linearities of up to 
about 10% there. 

Because our analysis of thermopower makes extensive use of p data, examples of the 
temperature dependence of p on both the metallic and insulating side are shown in Figs. 1 
and 2. They are very similar to those seen in previous work.00 

Well into the insulating regime (roughly n < 0.65 x lO^^m^^) p can be well represented 
by an equation developed for variable range hopping (VRH) through localized states in the 
presence of a Coulomb gap at the Fermi energy Ep (referred to as Efros-Shklovskii VRH), 

p = p,exp(TWT)i/2, (1) 

where pc and Tes are constants that depend on n but not T. Examples are given in Fig. 1 
for Sample 1. The values that we obtain for the constants are similar with those determined 
by others.lll Tes varies from ~ 35 K at n = 0.650 x lO^^m^^ to ~ 90 K at 0.48 x lO^^m'^. If 
the exponent is left as a variable we find an average value of 0.50 ±0.05 over this same range 
of n. We find poorer fits assuming Mott VRH (no Coulomb gap at Ep) which corresponds 
to a power law of 1/3 in the above expression. 

For n > 0.70 x lO^^m"^, a simple activated behaviour 

p = PcBxp (Ta/T) (2) 

provides a very good fit to the data below about 1 K. Thus these results suggest a transition 
from simple activated behaviour to VRH as n decreases. 

In the metallic region {n > uq) we have fitted our data to the equationlli 

p = Po + piexp[-{To/Ty] (3) 

where po, pi, Tq and q are constants that depend on n but not T. Examples are shown in 
Fig. 2 for Sample 2. For n > 1.5 x lO^^m"^, excellent fits were obtained by fixing q = 0.5. 
For lower n, the values of q (if this is left as a variable) increased and approached 2.0 at 
n 1.0 X lO^^m"^. The model curves fit the data successfully over the range 0.3 K (our 
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lowest T) to about 0.3 Tp. {Tp = 7.3nK with n in 10^^ m ^ units). Above this region the 



since Tq decreases as uq is approached (Fig. 2 indicates Tq for the lowest densities), a lowering 
of the upper temperature limit is also consistent with this. The reduced range of fitting at 
low n means that, pi, Tq and q are not well known here but, even so, po retains a relatively 
low error. 

Figure 3 shows the results on ctq = 1/po as a function of n.0 For all three samples (Tq 
follow the critical behaviour 



The solid lines are the best fits with the following parameters, with a in units of e'^/h (the 
data points at the nominal values of hq are (t(0.3K) and were not used in the fitting). For 
sample 2, a-m = 0.0 ± 0.3, as = 13.6 ± 0.7 and u = 0.94 ± 0.03. A good fit was not possible 
over the whole range of n and so was restricted to n < 3.0 x 10^^ m~^ (a range of [n/riQ — 1) 
of about 0.03 to 3). Sample 1 has an identical behaviour within experimental error. The 
higher-mobility Sample 3 also follows the same equation with am = 0.36 ± 0.15, = 34 ± 5 
and u = 1.39 ± 0.05. We have a smaller range of n for this sample, but the fit still covers 
a range {n/riQ — 1) 0.02 — 0.4. These results suggest u increases with peak mobility but 
clearly more data on a variety of samples are required. The values of a^, for Sample 1 and 2 
are consistent with zero within experimental uncertainty. The fits are good over the whole 
range except for the last points near n = hq which do not fall on the curves. This is not 
necessarily unexpected because a value of a = is not possible at a finite temperature due 
to thermal broadening of the electronic distribution. Small errors in no will have a similar 
effect. For Sample 3, a^, may be finite. However, if no is allowed to decrease from 0.96 
to about 0.93 x 10^^ m~^, a fit which is essentially indistinguishable over the range of the 
data can also be obtained with am = 0.0 ± 0.2, o"s = 32 ± 5 and u = 1.48 ± 0.05; both 
fits are actually shown on Fig. 3. A small discrepancy in uq could easily arise from the 
identification of the critical density for the MIT with that density, no, where dp/dT changes 



curves usually exhibit a maximum.lillil 



To also provides a natural scale for the data and. 




(4) 
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sign, a procedure which has no firm physical foundation;lli in addition there are uncertainties 
and possibly non-linearities in the determination of n from gate voltage values as mentioned 
above.0 

The critical behaviour described by Eq. (^) with am = is formally the same as that 
expected for a (continuous) Anderson transition with a 'mobility edge' at hq, whereas a 
finite am would correspond to a (discontinuous) Mott-Anderson transition; neither transition 
should arise in a non-interacting 2D gas.S'i Our p data in the insulating region are also 
consistent with this scenario in the sense that we observe activated behaviour for n below 
uq. Similar critical behaviour, usually with am consistent with zero, has been seen in many 
3D systems, typically with valued! of u in the range 0.5 — 1.3. There are also two previously 
reported cases related to 2D. Hanein et a/.@ have made a similar analysis to the one above 
for a 2D hole gas in GaAs and found a linear relation between ctq and n, but with a finite 
am- Feng et alM also found critical behaviour in ao in a low mobility Si-MOSFET though 
we note that their data all show negative dp{T)/dT, even in the range which they identify 
with 'metallic' behaviour. They found a{T) increased with T and fitted a dependence. 
This behaviour was attributed to local momemts and is so different from that found here 
that the two cases may not be closely related. Recently^ a model based on percolation of 
non-interacting electrons through local quantum point contacts also led to zero-temperature 
conductivity consistent with our observations. 

Note that the data presented in Fig. 3 were obtained from extrapolating a to T = 
by a procedure which focusses only on the 'strong' exponential, T— dependence of Eq. (3). 
It ignores any other weaker dependences^ that may be present below T = 0.3K, including 
those due to weak localization and screening. 

We now turn to the thermopower data. A selection of data on 5 as a function of T 
is shown in Fig. 4. From our previous work0 we know that the diffusion thermopower, 
5"^, is almost zero at n = 8.5 x 10^^ m^^ and one sees only phonon drag, which varies 
approximately as at the lowest temperatures. The fact that S'^ is very low at high n 
predicted by Karavolas and ButcheiH who attributed it to a particular combination of 
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scattering mechanisms. Indeed, when n 3> no the behaviour of both S"^ and 5*^ appears to 
be in generally good agreement with theory.0 

The dependence of obtains only in the Bloch limit, i.e., when q <^ 2kF where q is 
the average phonon wave number and kp the Fermi wave number. At n = 8.5 x 10^^ m~^ 
the condition q ^ 2kF is satisfied below about 1.4 K. Above this temperature one sees a 
gradual decrease in the exponent to roughly oc by 4K, but no simple power law is 
expected in this region. 

As n decreases, S begins to show two distinct regions with different T dependences. 
At T < 1 K, 5* has a much weaker, approximately linear, T-dependence indicative of S'^ 
becoming dominant. For n < uq, this low-T dependence, which is characteristic of ordinary 
metals, is replaced by an upturn in S; we will return to this feature later. 

Concentrating first on the metallic region, we analysed the data over the range 0.3-1.5 K 
to obtain S'^ by assuming 

S = S'^ + S'^ = aT + pT' (5) 

where a (= S'^'/T) and /3 are constants that depend on n. For the 2 highest densities we 
obtained good straight lines by plotting S/T as a function of T*~^ with s = 6 as expected. 
At lower n, proportionately lower temperatures are required to reach the Bloch limit and 
so we no longer expect to see oc in our temperature range. Hence we simply used 
s as a variable to give the best fit. For n = 1.9 — 5 x 10^^ m~^ we found s = 4, and 
for n = 0.97 — 1.5 x 10^^ m~^, s = 5. That s seems to increase again at the lowest n is 
unexpected, but more data are required before this is taken to be an experimental fact. 
We have no reason to believe that is not well-behaved in this temperature, density and 
resistivity range, but we must caution that our understanding of in these circumstances 
is not well founded as yet .0^1 Given this it is important to note that the value of a depends 
only weakly on the choice of s. Indeed, providing has a much stronger T dependence 
than linear, then simply taking the measured values of S/T at 0.3 K to be equal to a gives 
values in good agreement with those obtained using Eq. (|^). 
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Fig. 5 shows that a as a function of n appears to diverge as n ^ hq. One would expect 
such a divergence if Ep approaches a band edge or a gap in the DOS giving 5'°' oc l/ra. Das 
Sarma et al^ have suggested that a similar mechanism, carrier freezeout, is responsible for 
the MIT in 2DEGs. However, the present results are inconsistent with these explanations 
because Hall data, both our own and those in Ref. |TB], show that in the vicinity of the 
mobile carrier density is not approaching zero but equals n within a few percent; also any 
temperature dependence of n at fixed Vg is very small in the range 0.3 — 4.2 K. Given that 
theory successfully predictedil S'^ to be zero at high n as mentioned above, it is possible 
that the present results have an explanation in terms of the specific scattering processes 
important at low n, but there is no theory available for this situation at present. 

However, Eq. also implies a divergence of S'^. With the assumption of a constant 
density of states (DOS), Eq. (P is consistent with 

aiEp) = cr,n + (^s iy-^ - 1 j • (6) 

Again, with am = this is formally equivalent to an Anderson transition with Ec being the 
mobility edge. Taking am = 0, the use of the Mott relation 5"^ = - {n"^ k'^T /3e) {din a /dE)Ep 
with Eq. (|) then giveJl^li 

3e{EF-E,)- ^'^ 

The use of the Mott relation, and hence this result, is valid only if [Ep — Ec) / ksT ^ 1; in 
the opposite limit S'^ tends to a constant0ii (~ 228 fiV/K in 3D). The saturation is a direct 
consequence of the fact that the contribution to the conductivity of electrons below E^ is zero 
in this model. The relevant integral00 for S'^ has a weighting factor a{E){E — Ep) which, 
under normal conditions, leads to 5'"' being a measure of the derivative of a{E) with respect 
to energy, and hence to the Mott relation. This is consistent with Eq. (|^). However, when 
Ep = Ec only electrons above Ec contribute to the integral, which leads to a constant value 
of S'^. Numerical calculationslil show that the approximation of Eq. ([7p gives a magnitude 
roughly a factor of 2 too large when {Ep — Ec)/kBT ^ 2 (i.e. because of the approach to 
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saturation). For our sample at T = 0.4 K this corresponds to {n — no)/no = A (say) ~ 0.11, 
using the ideal DOS, go, with an effective mass of 0.19 mo. If the present model was indeed 
appropriate, it suggests that we should see saturation of S'^ as T increases for the samples 
with n = 0.97 and 1.06 x 10^^ m~^ (for higher n any saturation would be masked by the 
rapidly increasing contribution from S^). We are unable to see any obvious evidence for 
saturation at these densities and Eq. (|^) still seems to be valid there. However, we should 
recall that is not yet understood^ at low n and high p and there remains the possibility 
that an unexpected behaviour of could mask important features in S"^. 

Proceeding with the model, we have fitted Eq. (0) to the measured a, allowing for the 
possibility of saturation as hq is approached. To simulate a saturation we add A in the 
denominator (but allow it to be a variable when determining the best fit to the data) and, 
rewriting Eq. (^ in terms of n, we have 

a = S'^/T = -K/Ja^ + (— - 1)2 (8) 
V no 

where K = vn'^k'^/ {"ieEc) is a constant expected to be about 32yuV/K2 for Sample 1, again 
using go- If am is finite in Eq. (|]), then the Mott relation suggests that it will contribute 
to the denominator of Eq. (|^), also softening the divergence at n = no- However, the 
experimental am is so small that this is probably negligible compared to the finite-T effect 
considered here. 

The best fit of the data to Eq. (P) is shown by the solid line in Fig. 5 and uses A = 
0.15 ±0.01 andi^'= (9.5 ± 1.5) yuV/K^. Only data with n < 4 x 10^^ m'^ are fitted, a range 
consistent with that used for ao- A is consistent with that expected from the argument 
above, but K is too small by a factor of about 3. As with a, the fit can be improved if no 
is decreased, perhaps by 0.05 x 10^^ m"^, in which case A decreases and may even become 
zero. 

In spite of the relatively good correspondence between the model and the data, we em- 
phasize that we should be cautious in necessarily concluding that the model is fundamentally 
correct. We are comparing our results on a 2D system with a theoretical model of an Ander- 
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son MIT believed to be appropriate only to non-interacting electrons in 3D. At this point it 
is difficult to know whether the situation in 2D might be significantly changed by including 
interactions along with the strong disorder; this is a complex and ongoing theoretical prob- 
lem (e.g. see Refs. 0,^and references therein). Some progress has been made on calculating 
S'^ with the inclusion of weak interactions and disorder.!! We are, however, unable to ex- 
plain the observed strong n dependence in terms of the calculations which predict universal 
values. The calculated corrections are logarithmic in T. These would be the equivalent of 
similar corrections in p and, if they exist, would require much lower temperatures than we 
have available and would be difficult to detect. On the other hand, these facts mean that the 
analysis of a, po, and their relationship as zero-temperature quantities, is a self-consistent 
procedure which has a semi-classical meaning. Nevertheless, whether or not a goes sharply 
to zero near hq as would be required in the Anderson model, the observed change in be- 
haviour of S'^ as n is varied around no points to a fundamental change in transport taking 
place in this region. 

We should finally mention that we can also represent the data reasonably well over the 
same range using the simple expression a = — /iV/K^, with n in units of 10^^ m~^, 
but this has no obvious physical explanation; in particular, it does not have the form that 
we would expect for S'^ approaching a band edge at n = 0, i.e. a oc 1/n. 

Our data in the insulating regime also show a behaviour qualitatively consistent with 
activation across an energy or mobility gap. In such cases we expect an upturn of S'^ at low 
temperatures according to 

where A is a constant of order unity. In the particular case of the Anderson modeli'Bii 
Aks/e corresponds to the saturation value when Ec = Ep noted above in our discussion of 
the metallic region. When n is close to, but less than, no, we see activated behaviour in 
p (see Eq. (|^)) which is in accord with this. We do not see a high temperature saturation 
of 5* as Eq. (||) predicts, but this is not surprisingly because it will be masked by other 
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contributions to S, particularly 5"^ (see also below). 

There might also be another contribution to from VRH. Demishev et a/.ll have re- 
cently demonstrated the existence of such a component in S"^ using a 3D GaSb sample. The T 
dependence of our p data and other previously published data0 well into the insulating region 
are consistent with Efros-Shklovskii VRH. (When two or more conduction mechanisms are 
present, the appropriate S'^ are weighted by their contributions to cr). For this mechanism 
one expects S"^ to be a constant given bylil§l S"^ = -{kB/e){kBTEs/C){d\ng{E)/dE)Ep 
where Tes is obtained from the temperature dependence of p in Eq. (|lD, g{E) is the back- 
ground DOS, i.e, on which the Coulomb gap is superimposed, and C a dimensionless constant 
^ 6. If we take {dhi g{E) / BE) Ep ~ IZ-E'f (implying that Ee may be in the tail of the DOS) 
and again using to estimate Ee-, we find that the calculated S*^ are typically a factor of 
two smaller than the values of S at the observed minima. The argument is not significantly 
changed if Mott VRH is assumed.iii§ In this case S'^ oc T^/^ but the magnitudes calculated 
for S'^ are similar. The references give the relevant details. 

In the insulating regime, one expects a diverging localization or correlation length as the 
critical point is approached. Given this, one might also question whether the thermopower 
results are in the linear region close to the critical point. However, on this side of the 
transition we are never very close to the critical density. The closest point is at about 
uq — n 0.5 X lO^^m"^. The correlation length estimated in Ref. S under these conditions 
is about 20 /im. Assuming a magnitude of the disorder potential of order 1 K gives a threshold 
electric field for non-linear effects of order 80mV/m whereas the thermoelectric electric field 
is about ImV/m over the range of 0.7-7 K. Thus non-linear effects do not appear to be a 
problem for the present data. 

To put the present results in perspective, as far as we are aware the only previous work 
which attempted to follow 5"^ into the region of 2D electron localization was that of Burns 
and Chaikinil on thin films of Pd and PdAu. They found an upturn of 5'°' in the strong 
localization region but no divergence at higher conductivities. The authors attributed their 
results to the opening of a Mott-Hubbard gap. In 3D, Lauinger and Baumannil observed 
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critical behaviour of cr and a divergence of S'^ for metallic AuSb films, but the magnitude 
of the latter was 2 orders of magnitude smaller than seen here. Other experiments on bulk 
Sipi and NbSiii saw no divergence of S'^ on the metallic side. 

We close this section by making a few comments about S at higher T. appears to 
be dominant in this region because the strong T dependence of S (roughly to T^) is 
inconsistent with any other mechanism. As we have already indicated, little is known about 
the behaviour of where n is low and p is high,lliii and these are just the conditions that 
pertain around no. Thus, although should be present on the metallic side, its precise 
form is not known with any certainty. 

When conduction is dominated by VRH, we expect = because drag is based on the 
conservation of crystal momentum for electron-phonon scattering which will not hold for 
transitions between localized electron states.S'H Thus, on the insulator side should only 
exist when carrier excitation to delocalized states occurs, though we emphasize that there 
are no calculations appropriate to these conditions. Our p data are consistent with activated 
behaviour for n just below uq, but deeper into the insulating state VRH becomes dominant 
and so we would expect to see a strong diminution of in this region. Figure 6 shows 
experimental data on S* as a function of n at fixed temperatures of 2 K, 3 K and 4 K. We see 
that S rises as n decreases but crosses no smoothly, i.e., we no longer see divergent behaviour 
of S at riQ as we did for S'^. Indeed there is no feature that indicates that anything unusual 
occurs at hq. This behaviour is consistent with activated conduction below uq. However, we 
continue to see a rise in S, and presumably S^, to the lowest densities, which is not expected 
from our argument above. This can only be understood if significant activated conduction 
is also present at all densities. In this light, the fact that our p data at lowest n in Fig. 2, 




appear to follow the Efros-Shklovskii VRH model so well over a 



limited range of temperatures appears to be somewhat coincidental. 
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IV. SUMMARY AND CONCLUSIONS 



Generally speaking, the main focus of attention with regard to the MIT in 2D systems 
has been whether the metallic behaviour is the result of a transition to a new state induced 
by the strong disorder and possibly electron-electron interactions, or whether it is the result 
of a conventional physical mechanism that has yet to be unambiguously identified. 

On the metallic side, our analysis has mainly concentrated on data at low T, essentially 
below the region of rapid T dependence of p. Our main result is that the critical behaviour 
of ctq and S'^ of Si-MOSFETs in this region are surprisingly consistent with equations that 
are formally equivalent to those describing a 3D Anderson MIT. On the insulating side, the 
behaviour of p and S"^ over a wider range of T are also consistent with such a scenario. 
However, not all features that we see can be understood in this way. In particular, although 
the Anderson model exhibits scaling behaviour in both cr(T) and S'^{T)^ it does not appear 
to be able to explain the large increase of p with T in the metallic region observed in Si- 
MOSFETs. 

An Anderson transition can arise purely from disorder in 3D but general scaling results 
predict that the equivalent transition should not occur in 2D. Still, given that many features 
of the data do mimic an Anderson transition, a key question in the present work would then 
be whether a 'mobility edge' actually exists, in particular at low or zero temperature. If we 
look for more conventional explanations of the data, a point to bear in mind is that kpl ~ 0.3 
at no, where / is the electron mean free path estimated from the conductivity, which implies 
very strong disorder. Under these conditions it seems plausible that a sufficiently rapid drop 
in p with decreasing n could be mistaken for a mobility edge, but it might actually be caused 
by a smooth, though very rapid, transition from weak to strong (exponential) localization 
of the carriers in keeping with the scaling model. If this is the explanation, then it remains 
to be shown how this rapid drop can mimic a critical behaviour over a relatively wide range 
of n. 

In our analysis, the apparent divergence of S'^ as no is approached from above is basically 
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a reflection of the rapid drop in ctq, a result that might remain vahd even if no Anderson 
transition occurred. An unambiguous indicator of a mobihty edge would be saturation of 
S'^ as a function of temperature. Our data near no show no visible indication of this, but 
the presence of phonon drag is a complicating factor here and might mask any saturation. 
Clearly it would be an advantage to suppress phonon drag and reveal S''' over a wider range 
of T, but this is difficult to do in 2D. In principle, it could be done by the use of very 
thin substrates. We do not see any features in the overall behaviour of S which appear to 
correlate with the strong T dependence of p in the metaUic region. 

Finally, an important point that must be stressed is that to reliably identify any observed 
critical behaviour with a MIT requires data in the zero T hmit. Although our analysis of 
(To and S''' is based on an extrapolation to zero T, the actual data extend only to 0.3K. 
If the strong drop of p is caused by conventional physics, then we are examining the low 
temperature limit of this mechanism, but not necessarily the true low temperature limit of 
the system. This is true of practically all the experimental data published so far. 

In spite of these reservations, it is clear that the observed critical behaviour in both ctq 
and S''' indicates an unknown but interesting physics of strongly disordered and interacting 
systems. The results open up a new window on these systems that further constrains any 
theoretical model proposed to explain the MIT, whether such a transition be a quantum 
property or simply the result of a conventional physical mechanism. 
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FIGURES 

FIG. 1. Examples of resistivity p at various densities (in units of lO^^m"^) as a function of 
temperature T for Sample 1 in the insulating region. The data are plotted in the form log p as a 
function of T^^/^ to show that Eq. (|^) (Efros-Shklovskii VRH) gives a good representation at low 
n (0.52 and 0.61 on the figure). At higher n (0.70, 0.79 and 0.88) the data are well described by 
Eq. (^) below 1 K, corresponding to simple activated behaviour. 

FIG. 2. Resistivity p at various fixed densities n (in units of lO^^m"^) as a function of tem- 
perature T for Sample 2 in the metallic region. The solid lines are fits to Eq. (^) using only the 
points represented by open symbols. The closed symbols for the data near no are points not used 
in the fitting (see text). The short vertical lines give the value of Tq for each density. The data at 
n = 0.96 are nominally at no. 

FIG. 3. The main panel shows the density dependence of the conductivity ctq in the T ^ 
limit of Sample 2. The inset shows the same data for Sample 3. In both cases circles correspond 
to data obtained using Eq. (^). The lowest points designated by triangles are simply the measured 
cr(T=0.3K). The dashed vertical lines are the values of no used in obtaining the fitted curves. In 
the case of Sample 3, two fitted curves are shown (but are almost indistiguinshable over the data 
range) corresponding to the two values of no. 

FIG. 4. The thermopower S for Sample 1 at various fixed electron densities n (in units of 
-|^q15 j^-2j fuuctiou of temperature. 

FIG. 5. Density dependence of a = S'^/T for Sample 1. The closed symbols are obtained 
using Eq. (^) and the open symbols are simply the measured S/T at at 0.3 K (the values of a from 
the two methods at n > 3 x 10^^ m~^ are indistinguishable). The points at n = 0.97 x 10^^ m~^ 
are just below no but ksT broadening should make these indistinguishable from ng. The line is 
the best fit to Eq. (§) for n < 4 x 10^^ m.-^. 

FIG. 6. Thermopower at fixed temperatures of 2.0 K, 3.0 K and 4.0 K as a function of density. 
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